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Abstract 

A class of exact solutions of the dispersionless Toda hierarchy con- 
strained by a string equation is obtained. These solutions represent 
deformations of analytic curves with a finite number of nonzero har- 
monic moments. The corresponding r-functions are determined and 
the emergence of cusps is studied. 
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1 Introduction 



Integrable contour dynamics governed by the dispersionless Toda (dToda) 
hierarchy is a multifaceted subject. It underlies problems of complex analysis 
interface dynamics (Laplacian growth) j3] , Quantum Hall effect |lj 
and associativity (WDVV) equations 5\. A common ingredient in many of 
its applications is the presence of random models of normal N x N matrices 
H El 5 13 with partition functions of the form 

Z N = J dMdM t exp(-itr W{M,M^)), (1) 

where 

W(z, z) = zz + v - J2(tk z h + h z h ). (2) 
fe>i 

In an appropriate large N limit (H — > 0, s :— hN fixed), the eigenvalues of 
the matrices are distributed within a planar domain (support of eigenvalues) 
with sharp edges, which depends on the parameters t := (s = s, ti,t 2 . . .). 

If the support of eigenvalues is a simply-connected bounded domain with 
boundary given by an analytic curve 7 (z — z(p), \p\ = 1), then (s,ti,t2 ■ ■ •) 
are harmonic moments of 7 and the curve evolves with (t, t) according to 
the dToda hierarchy. Moreover, the corresponding r-function represents the 
quasiclassical limit of the partition function A particularly interesting 
feature is that for almost all initial conditions the evolution of 7 leads to crit- 
ical configurations in which cusp-like singularities develop. This behaviour 
is well-known in Laplacian growth [Bj and random matrix theory [Hj. 

In order to obtain solutions of the dToda hierarchy describing contour 
dynamics one must impose a string equation which leads to a particular type 
of Riemann-Hilbert problem ^ n ^ ms paper we present a method for 

finding solutions in the form of Laurent polynomials 

1 . U K-1 / Q \ 

z = rp + u + ■ ■ ■ + -jrzr, (3) 

which describe dynamics of curves with a finite number of nonzero harmonic 
moments, namely tk = tk = for k > K. We exhibit examples for arbitrary 
K and derive their corresponding r-functions. Furthermore the emergence 
of cusps is analytically studied. 
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2 dToda contour dynamics 



Let z = z{p) be an invertible conformal map of the exterior of the unit circle 
to the exterior of a simply connected domain bounded by a simple analytic 
curve 7 of the form 

z = S(z), (4) 

where bar stands for complex conjugation (zip) = z{p~ l ) on 7) and the 
Schwarz function S(z) is analytic in some domain containing 7. 
The map z(p) can be represented by a Laurent series 



00 

Z (p)=rp + J2^, (5) 
k=0 y 

with a real coefficient r. The coefficients (r, Uq, Ui, . . .) are functions of the 
harmonic moments t = (s = s, ti,t 2 . . .) of the exterior of 7, which in turn 
can be introduced through the expansion of the Schwarz function 

^HE^^ + i + E^Ti' ( 6 ) 

k=l z k=l z 

with (vi, V2, ■ ■ ■) being functions dependent on t. As a consequence of (HJ-©, 
it follows that z(p, t, t) solves the dToda hierarchy 

d tk z = {H k , z}, di k z = -{H k , z}, 

(7) 



1 ... 1 . 



H k : = (^) >1 + ^(z fc ) , := (z fe )<_i + Uz k 



where {/, g} := p (d p f d s g — d p g d s f), the function z(p x ) is defined by the 
Laurent series 

00 

z(p- 1 ):=~ + J2ukP k , (8) 



* k=0 



and the symbols (. . .)>i ((. . -)<-i) an d (• • -)o mean truncated Laurent se- 
ries with only positive (negative) terms and the constant term, respectively. 
Furthermore, this solution satisfies the string equation 

{zip), zip- 1 )} = I. (9) 
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These properties can be proved through the twistor scheme of Takasaki- 
Takebe [3]. It uses Orlov-Schulman functions of the dToda hierarchy 

OO CO 

k=l k=l 

(10) 

co _ 

k=l k=l 

and can be summarized as follows: 

Theorem If (z, m, z, m) are functions of (p, t, t) which admit expansions of 
the form do]) . (jHj) . flTU|) and satisfy the equations 

m 

z = — , m = m, (11) 

then (z, z) is a solution of the dToda hierarchy constrained by the string 
equation Q. 




3 Solutions 

Equations (fTTjl are meaningful only when they are interpreted as a suitable 
Riemann-Hilbert problem on the complex plane of the variable p. Thus (z, m) 
must be analytic functions in a neighborhood D = {|p| > r} of p = oo and 
(z, fn) must be analytic functions in a neighborhood D' = {\p\ < r'} of p = 0. 
The statement of the Theorem holds provided A := D f] D' ^ 0. 

We next prove that equations (fTTj) have solutions satisfying and 
ffTUf) with 

In this way we assume 

m 



m 



0, k > K, t K ^ 0. 



K 



nt k z k + s + Y 



Vk 



k=l 
K 



k=l 



(12) 



k=l 



k=l 
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Given two integers r x < r 2 we denote by V[r!,r 2 ] the set of Laurent 
polynomials of the form 

c ri p ri + c ri+ i p ri+1 H h c r2 p r2 . 

Let us look for solutions of (jllj) such that z and z are meromorphic 
functions of p with possible poles at p = and p = oo only Then, as a 
consequence of the assumptions ©,© and (fT2j) . from (fTTj) it follows that 

zeV[l-K,l], zeV[-l,K-l]. (13) 

The equation m = m is equivalent to the system: 

m>i = m>i, (14) 



m = mo, (15) 
m<_i = m<_ x . (16) 



If we now set 



^nt fc (* fc )>i + mo + A; (17) 

with 



m = m 

fe=l k=l 



K 

m = s + ^ki k (z k ) . 
k=l 

it can be easily seen that m has the required expansion of the form (|12|) 
provided z and z satisfy (0) and (jSJ) . On the other hand, the expression (|17|) 

777- 

for m has an expansion of the form (JT2jl if the residue of — corresponding 

z 

to its Laurent expansion in powers of z verifies 

771 

ites(— ,*) = s. (18) 

Hence the problem reduces to finding z and z satisfying (|5)i.(|8)i.(|I8)i and 

m 

,= 7 . (19) 
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In view of (|13|) we look for z and z of the form 

z = rp + u H h —f^—r-, 

(20) 

Z = - + U ~\ h U K _ip K ~ l . 

V 

Now, in order to prevent z from having poles different from p 
we have to impose 

m(pi) = 0, 

where Pi denote the K zeros of 

r + u p-\ h u K ^i p K = 0. 

In this way by using the expression (|17|) of m, the only variables appearing 
in (J 19)) are 

(P, t, t, r,u ,..., w , ... , %_i), Wi := Ui. 

Thus, by identifying coefficients of the powers p\ i = 1—K, . . . , 1 we get K+l 
equations which together with the K equations ()2*T|) determine the 2K + 
1 unknowns variables (r, u , . . . , Uk-i, wo, ■ ■ ■ , wk-i) as functions of (t,t). 
Moreover, provided r is a real coefficient, the equations (fTTj) are invariant 
under the transformation 

T/(rt = /(i). 

Hence if (r, u , . . . , Uk-i, Wq, ■ ■ ■ , Wk~i) solves (|T9|) so does 

(r, w ,..., w K -i, u ,..., Sjf-i). 

Therefore, if both solutions are close enough, they coincide and consequently 
Wi = Ui, as required. 

To complete our proof we must show that (|TSj) is satisfied too. To do that 
let us take two circles 7 (\p\ = r) and 7' (|p| = r') in the complex p-plane and 
denote by T and V their images under the maps z = z{p) and z = z{l/p), 
respectively. Notice that due to © and (jHJ), the curves Y and V have positive 



and p = 00 
(21) 
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orientation if 7 and 7' have positive and negative orientation, respectively. 
Then we have 

Tfl Tfl 1 / Tfl I / Tfl 

Res( — , z) — Res(— , z) = (p — dz <b — dz 

z z 2%ix J v z 2m J T i z 



= <£ zd v z dp (p zd„z dp = <£ dJzz) dp = 0, 

2 ? 7T/ 7 P 1 2m Jy P V 2i7T/ 7 PV ' 1 

where we have taken into account that the integrands are analytic functions 
of p in C — {0} and that 7 and the opposite curve of 7' are homotopic with 
respect to C — {0}. Therefore, as we have already proved that m has an 
expansion of the form (fT2"j) . we deduce 

m ,m 
Res( — , z) = Res( — , z) = s, 

z z 

so that (fTHj) follows. 

Let us illustrate the method with the case K = 2. The polynomial pz 
has two zeros at the points 



Pi = V= ' P2 = V= ' 

2u\ 2u\ 
and from (|2Tj) we get two equations which lead to 
— 2r t 2 w + 4r *2 Wo Mi + t t\ u Ui + 4r t 2 w m U\ — r iiM x — 4r t 2 M , u 1 

— s Mo u 1 — ti Uq Ui — 2 t 2 m "^l + r *i u i = 0' 

(22) 

—2r t 2 u + 2r t 2 U\+r t 1 w Wi + 4r t 2 u u Ui — r s ^ — r t 1 u u 1 

(23) 



-2rt 2 ul ul-2r 2 t 2 ul = 
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Identification of the powers of p in (|19|) implies 

p : —2r 2 t 2 + rui = 0, 



p : 2r t 2 uo — r t\ U\ — 4r t 2 uq U\ + Mo u x = 0, 

p : —2rt 2 u + 2r t 2 iti + r *i uq U\ + 4r t 2 uq Uq U\ — s u l — t\ Uq m 1 
-2 1 2 «o ^i - 4 r *2 + «i «i = 0. 

(24) 

Then by solving equations (j22J)-(J2I|) we get the solution: 



pVi 2^~st 2 _ ti+^tih (25) 

which corresponds to the conformal map describing an ellipse growing from 
a circle [6J 

Solutions for K > 3 

Exact solutions associated to arbitrary values of K can be found from the 
previous scheme. However in order to avoid complicated expressions , we set 

*1 — *2 = ' ' ' = *KT-1 = *1 = *2 = ' ' ' = *if-l = 0. 

and look for particular solutions satisfying 

Ui = U 2 = ■ ■ ■ = U K _ 2 = Ui= U 2 = ■ ■ ■ = UK-2 = 0, 

or equivalent ly 

z = rp+— ir -, z=- + u K - 1 p . (26) 
Under the previous assumptions and from (fT7|) we have that 



m = Kt K r K p K + 8 + K 2 i K r K ~ l UK-i + (27) 



Thus, we see that (|2*T|) leads us to a unique equation since from (J2*7|) it follows 
that m depends on p through . Furthermore, if pi satisfies z(pi) = 0, then 



K 



r 



Pi 



UK-l 
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Therefore, (|21j) becomes 

s - — — + (K — l)Kr K ~H K u K ^ = 0. (28) 

Uk-X 

On the other hand, it is easy to see that 

m Kr K t K / 2 K _ x su K „i-Kr K+ H K 

-p + K T t K + 



— 2 I K 1 ' 

- uk-i V / V 

consequently, by equating coefficients and taking (J2T)|) into account, we find 
that ()19|) leads to two equations only . More precisely 

Kr K t K 

p : r = — , 

uk-i 

(29) 



u. 



Then, we get three equations for the three unknowns r, Uk-i, which 

proves that there exits a solution of the form (J2HJ). In fact, by solving (pSjl- 
(12511 we find that 

Kt K r K - x , , 

2 = r P + — in — ' ( 30 ) 

with r satisfying the implicit equation 

K 2 (K - l)t K t K r^ K -^ -r 2 + s = 0. (31) 

Figures 1 and 2 show examples of the evolution of the curve z(p), \p\ = 1 
as s grows and tx is kept fixed. 



r-functions 

In P it was proved that there is a dToda r-function associated to each 
analytic curve z = z(p), \p\ = 1, given by 

2 logr = -\s 2 + sv Q - ^y^(t fc Vk + i k Vk), (32) 
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Figure 1: solution corresponding to K = 4 



Figure 2: solution corresponding to K = 10 

where Vk are the coefficients of the expansion (JBJ), and v is determined by 

dv o , 2 <91og r 

^ = logr, »b = -^-. (33) 
For the class of solutions (}30|) we have 

2 logr = -^s 2 + sv - ~(t K v K + i K v K ), (34) 
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and from and (|3U|) it follows that 



^- £z z K dz =^r- j z(p)z(p) K (r-(K- l)^-i)dp 

(35) 

(r 2 - s){Ks - (K — 2)r 2 ) 



2K(K- l)s 



On the other hand by differentiating ()34|) with respect to s and by taking 
into account one finds 

(i(-2)(i(-l)%|^-(^-2)rV^ 
t; -- S + S logr + 2((ir-l) 2 ^|^| 2 r^-H) 

(36) 

, (^-2)(s-r 2 ) / | (iT-2)r 4 



2{K-\)K \ (K -l) 2 K 2 \t K \ 2 r 2K - r 4 
Thus, (p31Jl - (j3T)]l and (J3T|) characterize the r-function of the curves determined 

by nam. 

Cusps 

The pictures of the curves associated with (|3T)J) - (f3*T]) show the presence of 
cusps at some value of s for each fixed value of tx ■ Indeed by using the 
parametric equation p = e i9 , (0 < 9 < 2n) for the unit circle, we have that 
cusps on the curve z = z{p) appear at points where z$ = 0, Zgg ^ and 
zggg/zgg has a nonzero imaginary part. Therefore a necessary condition for 

P — p(8) is 

Dz 

— (p) = 0, |p| = l. 

Thus from (|50|l we deduce 

p K = K(K-l)i K r K - 2 , 
which together with the condition \p\ = 1 requieres that 

r = (K(K-l)\t K \)~£*, (37) 
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at some value s = s(tjc)- But according to (|3T|) one finds that this happens 
at the value s given by 

8o = tj^(K(K-l)\t K \)-™, (38) 

which is the point at which the profile of both positive branches of r, as 
functions of s, develop an infinite slop (see figure 3). 




Figure 3: The positive branches of r(s) for K = 10 
Therefore, there are K cusps given by the roots 

which emerge when s reaches the extreme value sq of the domain of existence 
of the two positive branches of r as a function of s. 
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